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ABSTRACT 

We report first results from the Anglo-Australian Telescope Rocky Planet Search 

- an intensive, high-precision Doppler planet search targeting low-mass exoplanets in 
contiguous 48 night observing blocks. On this run we targeted 24 bright, nearby and 
intrinsically stable Sun-like stars selected from the Anglo-Australian Planet Search's 
main sample. These observations have already detected one low-mass planet reported 
elsewhere (HD 16417b), and here we reconfirm the detection of HD 4308b. Further, we 
have Monte-Carlo simulated the data from this run on a star-by-star basis to produce 
robust detection constraints. These simulations demonstrate clear differences in the 
exoplanet detectability functions from star to star due to differences in sampling, data 
quality and intrinsic stellar stability. They reinforce the importance of star-by-star 
simulation when interpreting the data from Doppler planet searches. The simulations 
indicate that for some of our target stars we are sensitive to close-orbiting planets as 
small as a few Earth masses. The two low-mass planets present in our 24 star sample 
indicate that the exoplanet minimum mass function at low masses is likely to be a flat 
a ~ -1 (for dN/dM oc M") and that between 15±10% (at a = -0.3) and 48±34% (at 
a = —1.3) of stars host planets with orbital periods of less than 16 days and minimum 
masses greater than 3 . 
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Introduction 



The planetary mass function is a key constraint which theories of star and planet formation 
must be able to accurately predict, if they are to be considered viable. While a consensus is 
emerging that planet formation, in general, is probably dominated by the growth of rocky cores via 
the accretion of dust particles and/or ice s (i.e. the "core accr etion" paradigm developed to explain 



the formation of the Solar System, e.g. iPollack et al.lll99d ). the priniary a lternative scenario of 



direct gravitational collapse remains a plausible mechanism (e.g. iBossI 120071 ) 



Three indirect lines of reasoning suggest that terrestrial-mass planets could orbit at least a few 
percent of stars. The first planet with a minimum mass of less than 10 - GJ 876d - was identified 
j ust 3 years ago. G J876d has a minimum (msmi) mass of 4.9 Mq and a probable mass of 7.5 
(jRivera et al.ll2005l ). In the years since, a further five Dopp l er exoplanets have been announced 
with minimum masses of less than 10 (jUdry et al.l 120071 : iMayor et al.l l2008l ) . And a further 



two low-mass exoplanets h ave been revealed by gravitational microlensing searches (jBeaulieu et al 
20061 : iBennett et al.ll2008l ). So there can be little doubt that exoplanets below 10 M® do indeed 



exist, though whether such planets are rocky ones (like the Earth) or icy ones (like Neptune and 
Uranus) is yet to be unambiguously d emonstrated. Second, gas-giant planets exist around at least 
6.5% of nearby main sequence stars (jMarcv et al.l l2005al ) and planet formation theories suggest 
that rocky cores and embryos should accompa ny such Jovian planets (jWetherill &: Steward Il989l : 
Goldreich et al.l 12004 : iKenvon &: Bromlevll2006l ). Finally, dusty protoplanetary disks are very com- 
mon around young T Ta.uri stars, which suggests the ubiquity of th e building blocks of rocky planets 
(jllillenbrand et al.|[l99a : lllaisch et aPboOll : iMamaiek et aklbood ). 



While not definitive, these observations are suggestive of the existence of rocky planets outside 
our Solar System. However, the theory of rocky planet formation remains less than solid. We 
do not know what fraction of stars form rocky planets, nor how many of those planets avoid 
dynamical ejection by the larger planets. Worse, we do not know if rocky planets with masses above 
1 quickly accr ete gas and volatiles t o become ice-gia r its lik e Neptune and Uranus. This has 
been predicted bv llda &: LinI (j2004l ) and iGoldreich et al.l (j2004l ). and would result in a planetary 
mass desert between 1 and 14 M^, as indeed we see in our Solar System. Theory does suggest 
that disks with higher surface-mass density than the minimum-mass solar nebula are expected 
to produce numerous rocky plai iets formed clos er to their parent star, which would make their 
Doppler detection more feasible (llda &: Linll2004l ). Moreover, planets of 1-15 may retain water 
in amounts th at are comparable to the silicates and iron-peak elements, resulting in a family of 
ocean-planets (jLeger et al.ll2004l ). 



What is clear, is that the detection of a statistically meaningful sample of planets between 5- 
25 M^, along with measurements of their masses, radii, and orbits, will be required to significantly 
inform these diverse theories, and drive forward our detailed understanding of how planets form. 
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2. The Anglo- Australian Planet Search 



The Anglo-Australian Planet Search (AAPS) is a long-term radial velocity program targeting 
the detection and parametrisation of exoplanets. The AAPS' main survey targets 250 southern Sun- 
like stars, and has been in operatio n since 1998. To dat e the survey has discovered 32 new exoplanets 



orbit i ng stars in the main sam p le (ITinnev et al 



2002 



Jones et alJl2002l . l2003al jbl. [2003; ICarter et al 



2001 



2002a 



200 



i 



2003 



McCarthv et alJl2004l : 



200,i 
3t al.l [ 



2006 



Butler et al. 12001 



O 'Toole et al.ll2007l : 



Bailey et al.l l2008l). AAPS' pre cision Doppler measurements are made with the UCLES echelle 



spectrometer (jPiego et al.lll99ll ). An iodine absorption cell provides wavelength calibration from 



5000 to 6200 A. The spectrograph point-spread function and wavelength calibratio n are derived 



from the iodine absorpt ion lines embedded on every pixel of the spectrum by the cell ([Valenti et al 



19951 : iButler et al.lll996l ) . This observing and analysis system has d emonstrated long term precisions 
of 3ms~^ for late-F, G, and early-K dwarfs brighter than V=7.5 ( Tinnev et al.ll2005l : iButler et al 



200l|) 



The UCLES dispersion through a 1" slit gives an effective resolution of A/AA=60000. A 
precision of lms~^ corresponds to displacements of 0.5 millipixel at our detector, and this, then, 
is the level below which we aim to hold the sum of all systematic uncertainties. Many subtle 
uncertainty issues enter at the level of 0.5 millipixel, such as determining the wavelength scale to 
9 significant digits, establishing the PSF shape at 1 part in 1000 (which varies dramatically with 
wavelength and seeing), measuring the charge transfer efficiency of the CCD, and calibrating the 
non-linear response of the quantum efficiency of the CCD. In addition, we specifically treat the 
Earth's telluric H2O and O2 lines that contaminate our spectra at levels below 0.5% intensity, and 
we eliminate all ghosts and defects from the spectrometer and CCD, on a pixel by pixel basis. We 
also employ a fully relativistic Doppler correction to the barycenter of the Solar System to remove 
the effect of the Earth's orbital motion and rotational spin. 



Our data proce ssing procedures follow those described bv lButler et al.l (jl996l . l200ll . l2006l ) and 
Tinnev et al.l (|2005l ). All data taken by the AAPS to date have been reprocessed through our 
continuously upgraded analysis system. The results presented in this paper arise from the current 
version of that pipeline. 



3. The Anglo- Australian Rocky Planet Search 

The detection of very low-mass exoplanets by the AAPS and other Doppler programs within the 
last 4-5 years has in large part been due to the dramatic improvements achieved in the measurement 
precisions these searches achieve. These have improved to such an extent that it is now clear that 
noise sources intrinsic to the parent star themselves are the limiting factor for very low-mass 
exoplanet detection. Characterisation of these noise sources (jitter, convective granulation and 
asteroseismological p-mode oscillations) has become an important focus of Doppler planet detection. 
A few obvious modifications to current observing strategies have emerged - (1) target low-mass 
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stars; (2) target chromospherically inactive and slowly rotating stars; (3) target high-gravity stars 
(where p-mode oscillations are minimised) and (4) extend the observations of stars over several 
p-mode fundamental periods, so that asteroseismological noise is averaged over. 

The Anglo- Australian Rocky Planet Search specifically seeks to focus on the last three of these 
observing strategies, in an effort to push to the lowest possible detection limits achievable with 
our system. That is, it focusses on the observation of bright stars known to be relatively stable 
to radial velocity jitter (e.g. Wright 2005). We ensure that all observations of each target are 
of sufficient length to average over the dominant asteroseismological period for ea ch target star. 



For this purpose, we employed the algorithms of lO'Toole. Tinnev &: Joned (|2008l ) to derive the 
maximum asteroseismological beat period (-Pmax) for each star, and then ensure that observations 
are at least this long. Observations of our Rocky Planet Search targets, which typically spanned 
15-20 minutes, easily achieved this goal. 

Finally, and most importantly, observations are carried out over long, contiguous observing 
runs of 48 nights. These long observing blocks are critical, as they allow us to integrate up Fourier 
power and suppress sidebands in the window function over many periods. The first of these Rocky 
Planet Search observing blocks was scheduled on the 3.9m Anglo-Australian Telescope (AAT) from 
2007 Jan 10 to 2007 Feb 26, spanning two full bright lunations and the dark and grey time in 
between. 

The power of this intensive observing approach is m ost simply and drama tically demonstrated 



by our detection of the short-period planet HD 16417b (jO'Toole et al.ll2008bl ) from 24 epochs ob- 
tained on that run, and the non-detection of the same planet in a similar quantity of similar quality 
data spread over a 2 year period. This bears out the simulations performed by other investigators - 
that to detect planets with amplitudes of l-to-3 times that of the net measurement uncertainty, it 



is ne cessary to acquire observations over many contiguous periods (see for example, iNarayan et al 



20051). 



3.1. Sample Selection 



The stars chosen for Rocky Planet Search observation on this first run were a subset of the 
AAPS main sample (j Jones et al.l l2003al ) that are bright (V<6.7), inactive (logR^j^< —4.8) and 
in the right ascension range 0-17^* with southerly declinations, suitable for observation in January 
and February at the Anglo-Australian Telescope. The sample all had primary asteroseismological 
period, Pmax, less than 800s. It includes 24 F, G and K stars. No metallicity criteria were used 
in this selection. Stars with known stellar companions within 2 arcsec are removed from the 
observing list as it is operationally difficult to get an uncontaminated spectrum of a star with a 
nearby con ipanion. One sta. r which satisfied th ese criteria (HD 75289) has a known exoplanet in a 
3.5 d orbit (jUdry et al.ll2000l : iButler et al.ll2006l ) of such high-amplitude (K=54ms~^) that it would 
have significantly complicated the possibility of detecting a low-mass planet. As a result it was not 
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included in the observed sample. Otherwise, stars were not rejected on the basis of having known 
exoplanets. HD4308 was selected for inclusion in this sample when the observing program was first 
proposed in 2005, and was retained for observa t ion ev en after a low-amplitude planet (K=4ms~^) 
was announced in a 15.1 d orbit bv lUdrv et al.l (|2006l ). 



The stars observed, and their properties, are summarised i n Table [H All of t hese stars have 
been included in large-scale studies of nearby solar type stars. iNordstrom et al.l (|200J) included 
them in a magnitude-lirn ited, kinematically unbiased study of 16682 nearby F and G dwarf stars. 
Valenti &: Fischeii (|2005l ) included them in a study of the stellar properties for 1040 F, G and 
K stars observed for the Anglo-Australian, Lick and Keck planet search programmes. Valenti &: 
Fischer used high signal-to-noise ratio echelle spectra (originally taken as radial velocity template 
spectra) an d spectral synthes i s to derive effective temperatures, surface gravities and metallici- 



ties whereas 



of Alonso et al 



Nordstrom et al.l (120041 ) used Stromgren photometry and the infrared flux calibration 



(|l996l ^. Both studies use Hipparcos parallaxes to obtain luminosities and make 
comparisons with different theoretical isochrones to derive stellar p arameters. To determi ne stel- 
lar masses and ages Va lenti & Fischer use Yon sei-Yale isochrones ( iDemaraue et al ll2004l ). while 



Nordstrom et al. (2004) use Padova isochrones ( Girardi et al. 2000 : [Salasnich et al. 200ol ). Both 



sets of derived parameters are consistent to within the uncertainties of these studies. 



3.2. Observations 

The observing strategy on each night was straightforward - to observe every star in the target 
list for at least 15 minutes (or as long as is required to obtain S/N~400 per spectral pixel given the 
available transparency and seeing conditions) on each of the 48 nights of observing. The number of 
observations actually obtained for each of the 24 targets is given in Table [2] and ranges from 18 to 
44, with a median of 31. The second column in this table provides the root-mean-square (RMS) for 
each target about the mean velocity, with four exceptions - the values given for the two stars with 
planets (HD 16417 and HD4308 - see below) are the RMS of the residuals after a best-fit planet 
was subtracted, and the values for the two components of the binary system a Gen (HD 128620 
and HD 128621) have had linear trends subtracted. The "Act. jitter" column sho ws the activity 



jitter derived using an updated version (Jason Wright, private communication) of the lWrightl (j2005l ) 



recipe, which is considered to pro duce jitter estimates good to ± 50%. The "Osc. jitter" column 



gives the oscillation jitter (from lO'Toole. Tinney &: JonesI l2008l ) for the average exposure time 



used for each target. To present these results graphically, histograms of the measured Doppler 
velocities are shown in Figure □ for each object (with the plots for HD4308, HD 16417, HD 128620 
& HD 128621 showing residuals as described above). Overplotted are Gaussians with full- width at 
half- maximum (FWHM) equal to the RMS (dashed line) and the activity jitter (dotted line) for 
each object. 

There is at least one case (HD 84117) where the observed velocity variability distribution 
appears to deviate significantly from a Gaussian distribution. We therefore plot the actual velocities 
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obtained for this object in Fig. [2l which show evidence for significant variabihty over the course of 
the 48n run, though with no obvious periodicity, nor with an obviously Keplerian shape. (No other 
star observed on this run shows a similar variability trend, indicating that the variability seen is 
not a systematic effect of our measurement system). Adding in additional AAPS data of similar 
quality taken since 2005 Jan 30 indicates that HD 84117 has shown excess velocity variability since 
2005 over that expected from activity jitter alone (with an RMS of 4.7 ms~^). The periodogram of 
HD 84117 shows essentially no power at periods of less than 40 d, though a complex, broad power 
distribution is seen at longer periods, i.e. longer than the time-span of the observations. We can 
therefore fit the data with no compelling single Keplerian. HD 84117 may either contain multiple 
planets, which will require intensive observation to disentangle, or may be a star with an unusual 
class of velocity variability. We can nonetheless say with confidence that it does not host a low-mass 
exoplanet in an orbit of less than ~30d. 

Thirteen of the 24 stars have observed RMS values consistent with their jitter estimates to 
within ±50% and all but two are consistent to within a factor of two - this is in line with expec- 
tations if the jitter estimates have a 1-a uncertainty of ±50%. The two outliers in this comparison 
are HD 1581 (F9.5V) and HD 23429 (KOV), where the observed velocity variation is much less than 
that predicted by the jitter model (by factors of 2.3 and 2.5, respectively), and HD 115617 (G5V) 
where the level of variability observed is larger than that predicted by the jitter model (by a factor 
of 2.2). The activity jitter estimate does not include velocity variability due to asteroseismological 
oscillations (which is in general much smaller than the activity jitter - see Table [2]) and more im- 
portantly, does not include the effects of convective granulation, which is (to date) unparametrised 
for any Doppler planet search stars. What can be concluded is that at very high velocity precisions, 
even "inactive" stars remain velocity variable at low levels. 



3.3. Planets in the Sample 



There are two objects in our target list that have detectable planets orbiting them: HD 4308 
and HD 16417. The latter was discovered as part of this observing campaign with P=17.24±0.01 d , 
e=0.20±0.09, and msini=22.1±2.OM0. It is presented in a separate paper (jO'Toole et al.ll2008bl ). 
to which we refer the reader for further details. 



HD 4308b was discovered by lUdry et al.l (|2006l ) and is a ~ 13 M® planet in a 15.6 d orbit around 
a G5 dwarf. It was observed 25 times during our campaign and the 15.6 d peri od is clearly evident 
in the 2D Keplerian Lomb-Scargle (2DKLS) periodogram (|Q 'Toole et al.l 120071 ) shown in Figure [3j 
The data from this run alone are indeed sufficient to a Keplerian fit which measures similar orbital 
paramet ers (see Fig. P to th e Udry et al. data. A scrambled- velocities false-alarm-probability test 
(FAP - iMarcv et al.ll2005al ) indicates that for 5000 trials, '^16% of scrambled velocities for this 



data set give a reduced chi-squared better than the solution. 



To refine the orbital parameters of HD 4308b, we have combined the velocities published in 
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Udrv et al.l (|2006l ) with our AAPS velocities - this includes both the observations obtained of 
HD 4308 in the Rocky Planet Search run of 2007 Jan-Feb, and other AAPS observations obtained 
with similar SNR between 2005 Oct 21 and 2007 Nov 27. These observations are listed in Table 
O The 2DKLS periodogram based on this combined dataset is shown in Figure [SJ The inset in 
this figure shows the corresponding slice through the 2DKLS at constant period, as a function of 
eccentricity, and demonstrates how eccentricity is not well constrained by this data set. 

In Figured] we show our best fit to the HD4308 data from the Rocky Planet Search run alone, 
while in Figure [6] we show the combination of all high-SNR data from the AAT and published 
HARPS data, and list the derive d parameters in Table H (The mass for HD4308 of 0.91±0.05 M© 
due to IValenti &: Fischeii (|2005l ) is adopted, as is a jitter estimate for HD4308 of 2.17 ms^^). 
We have determine a FAP for this planet of <0.02%, based again on 5000 scrambled velocity 
datasets. Though the eccentricity derived from the combined data set (e = 0.27 ± 0.12) is higher 
than that published by Udry et al., we n ote two points in rela,tion to this; first, th at it has been 



recently demonstrated by several studies (lO'Toole et al 



200S 



Shen Turned I2OO8I ) that there is 



a systematic bias against measuring zero eccentricities for low-amplitude planets; and second, that 
the uncertainty estimates for eccentricity produced by the least-squares fitting of Keplerians can 



serio usly underestimate the true eccentricity uncertainty represented by a data set (jO'Toole et al 



2003 ). Given this, we do not believe the differences between the Keplerian parameters for HD 4308b 
derived here and by Udry et al. are significant. 

Finally, we note that the residuals to the Rocky Planet Search data alone (Fig. S]) are suggestive 
of a further periodicity in that data set at 30-40 d. To examine the possibility of there being a second 
planet present in this data, we have constructed the 2DKLS for the AAT and HARPS velocities 
with the Keplerian of Fig. [6] removed (see Fig. [7|). The result is suggestive of power at period 
between 30-80 d. While potential planets at 32 d and 48 d periods can be fitted to this data, they do 
not do so with a significance that warrants a claim to have detected further planets in this system. 



4. Simulations 



The biases inherent in planet search observing and analysis strategies remain the largest single 
hurdle to a robust understanding of the formation processes that have built the 300-odd exoplanets 
known to date. Considerable work has been cond ucted on trying to eliminate these uiider ly- 
ing biases from Dopp l er exoplanetary statistics (e .g. 
Gumming et~aDl2008l . boosl . Il999l : lGumminel[20o3 ). 



O 'Toole et al.ll2009l : IWittenmver et al.ll200(il : 



To assess the selection functions delivered by the AAPS, we have been working towards a 



detai led understanding of Doppler noise sources intrinsic to stars (e.g. 
20081 ). We then use detailed object-by-object Monte-Garlo simulations (jO'Toole et al 



O 'Toole. Tinney Jones 
20091 ). to 



explore the biases introduced by these noise sources, when combined with our observational window 
functions. 
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The procedure we employ is to generate single Keplerians for a grid of periods, eccentricities 
and semi-amplitudes, which are then sampled at the actual observation time-stamps and which have 
noise added to them in line with the actual measurement uncertainties and stellar noise appropriate 
for each epoch. These simulated observations are then subjected to an automated detection process, 
enabling us to determine the range of period, planet mass and orbital eccentricity to which each 
data set is sensitive. The simulations are performed using the Keter and Swin burne supercompute r 
facilities. Details of the simulations and detectability criteria are described in lo 'Toole et al. ( 2009 ) 
- though with the modification to the procedure in that paper (which describes the analysis of 
purely artificial data), that the simulatio ns of actual data, used here, include added n oise terms 



due to intrinsic stellar velocity variability (jWright 



2005 



'Toole. Tinnev Jonesll2008l ^. 



To simulate the 48 night s of data obtai r ied on this first Rocky Planet Search run, we used a 
subset of the model grid from lO'Toole et al.l ([20091) with input periods on a logarithmic scale from 



log^o = 0.3,0.6,0.9,1.2 (i.e. periods of between 2d and 16 d); input eccentricities of 0.0, 0.1 & 
0.2; and planet masses of 0.01, 0.02, 0.05, 0.1, 0.2 and O.SMjup. The eccentricity range simulated 
was deliberatelylimited to the "near-circular" range of e=0.0-0.2, since (a) the majority of planets 
in such small orbits (i.e. <0.12AU) will have undergone significant tidal circularisation, and (b) 



the vast majo rity of detected Dop pler exoplanets with periods of less than 16 d have eccentricities, 

a). As our data is approximately evenly sampled (each star was 



e < 0.2 (e.g. iMarcv et al. 



20051 



observed on average once per night at approximately the same time), we have supplemented the 
above with additional simulations at the same eccentricities and masses, but with non-integer 
periods of 2.8, 4.3, 5.6, 8.7, 11.2, 17.4 and 22.4 days. Because the Keplerian function of a Doppler 
exoplanet is unable to determine inclination angle to the line of sight, i, no attempt has been made 
to deal with this in these simulations. As a result, whenever we refer to "mass" for a simulated 
planet in the following discussion, we are actually probing the msini " minimum mass" of that 
planet. For each (period, eccentricity,planet mass) point in this grid, 100 realisations are performed 
each with random noise added appropriate to the combined impact of the uncertainty at that time- 
stamp and the jitter and p-mode noise appropriate for that star (see Table [2|). Each star, then, is 
the subject of 21600 simulations. 

These simulations allow us to generate estimates of the detectability, for each star, of putative 
planets at each (mass, period, eccentricity) point in the simulations. In this context, we define 
"detectabili ty" at a given set of orbital parameters as the number of detections (using the detection 
criteria of lO'Toole et al.l |2009|) divided by the corresponding total number of realisations. We 
have then integrated this detectability over eccentricity since ( a) detectability is a pproximately 
constant with eccentricity at the values we have simulated (see lO 'Toole et al.ll2009l ). and (b) we 
wish to focus on our survey's sensitivity to planet mass and period. The result is a set of surfaces 
indicating detectability as a function of m sin i and period. These surfaces are shown in Fig. [H] for 
two examples (HD4308 and HD 53705), as well as for the whole sample of 24 objects. 

The specific examples shown for HD 4308 and HD 53705 display the general features seen in 
the simulations of all the targets: first that detectability is very high at periods longer than a 
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few days and at masses above 20 M®; and second, that the details of the detectabihty contours 
vary from object to object. The primary cause for this is the different time and SNR samphng of 
the observations of each star. This reinforces the importance of simulating radial-velocity data on 
a star-by-star basis, rather than on a survey sensitivity basis. We not that the structure in the 
surfaces at ~ 2 and ~ 4 days is an effect of sampling; the observations were taken on average once 
per night at approximately the same time. This means that planets with periods that are integer 
multiples of this have poorly constrained parameters and are therefore harder to detect. 

And finally, we must ask how many planets are detected from our actual observations if we 
apply to them the same automated criteria used to detect planets from the simulated data sets? 
Reassuringly, we find that just two planets are detected (HD 4308b and HD 16417b, although the 
former is perhaps marginal) - none of the other stars reveal planets that pass our automated 
detection criteria. 

5. Discussion 

We can see from the average detectabihty surface in Fig. [8] that the sensitivity of our survey 
of these 24 stars extends down to detectabilities of 10% or more at masses of 5 over periods of 
2-10 d. For higher mass planets, the survey is more sensitive, with the detectabihty being 40% or 
better for all planets with minimum mass > 14 and period > 3 d. 

With two planets known with minimum masses of 13.0 and 22.1 from a survey of 24 stars, 
the first-order conclusion that can be drawn from this data set is that for periods of 2-16 d the 
frequency of low-mass (i.e. msini=lO-25M0) planets is ~ 8 it 6%. However, this neglects both 
the fact that we know that our detectabihty is less than perfect (though non-zero) over this mass 
range, and the fact that we know that our planet detectabihty is a strong function of planet mass 
and a weak function of planet period. These effects need to be accounted for if we are to draw 
more robust conclusions from this survey. 

For different exoplanet mass functions, we would expect our survey to have sensitivity to dif- 
ferent masses. Put simply, a very steep exoplanet mass function should counteract our detectabihty 
surfaces to enable the detection of planets in the very lowest masses. Conversely, a fiatter mass 
function should bias our results toward the detection of higher-mass planets. So we should be 
able to use our data and simulations to provide much more insight into the frequency of low-mass 
planets. 

If we parameterise the planet mass function as a power-law dN/dA-I oc M" (we assume for sim- 
plicity that there is no change in the incidence of planets over this period range at these masses), we 
can ask, given our known detectabihty surface as a function of mass and period, what normalisation 
of that power-law would give us two exoplanet detections from our survey? Given our assumption 
that planet frequency is independent of period in this period range, and because detectabihty is 
also a weak function of period in this range, we make the further simplifying assumption that we 
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can integrate over the detectability surface at a given mass to derive an average detectability as a 
function of mass alone Dp{M). 

We can then write the number of planets detected from our survey, N^et, as a function of the 
number of stars observed N stars, the average frequency of planets in a given mass and period range, 
F, the detectability as a function of mass, Dp(M), and the power-law mass function. 



from which (given Np)et = 2) we can derive F (averaged over the period range P=2-22 d 
and the mass rangfQ 0.01-20 Mjup), which varies from 0.14±0.10 at a=-0.25, up to 0.70±0.85 
at a=— 1.75, in the F, G and K stars probed here. These planet frequencies might seem high in 
comparison to other work, though this i s almost certain l y due to the much lower masses being 



considered. For example, in Figure 12 of ICumming et al.l ([200a), those authors find an exoplanet 



frequency of 3% for periods less than 16 days and masses above 30 M^. The data presented here, 
in comparison, include significant sensitivity down to 4M0. The normalisations found here are, 
of course, uncertain due to the Poisson-counting statistics of just 2 detections (in a sample of 24 
stars), and do only probe a relatively small range of periods (or equivalently, semi-major axes). 
Having said which, this survey is important in being one of the first Doppler surveys of low-mass 
exoplanets with robustly characterised selection effects. 

Fig. [9] shows how the expected number of detected planets varies as a function of different 
underlying mass functions (along with the msini minimum masses for the two detected planets 
in this survey, HD 4308b and HD 16417b). As expected, a different peak in the expected number 
of observed exoplanets is predicted for different mass functions. a=— 1.75 leads to a peak in the 
number of detected planets around 0.02 Mjup (6.3 M®), the a=— 1.0 expected detections peak at 
around 0.05 Mjup (16 M®) while the flatter a=— 0.25 mass function produces an expected detection 
peak at around 0.10 Mjup (32 M®). 

Given that we have found two exoplanet signals in the 16-25 minimum mass range (cor- 
responding to the 0.05 Mjup or 16 bin in Fig. [9]) we may consider these as providing the first 
available rough limits on the underlying mass function as determined by how changes in the mass 
function cause the peak of detections to move toward and away from the 0.05 Mjup (16 Mq) bin. 
We find the peak moves from the 0.02 to the 0.05 Mjup bin for a values shallower than a=— 1.3 and 
from 0.05 Mjup to 0.1 Mjup at a=-0.3 (corresponding to F values of of 0.48±0.34 and 0.15±0.10 
respectively). Our detections thus indicate the mass function lies within these limits and that the 
exoplanet mass function at low masses is more likely to be roughly flat a ~ — 1 rather than steep. 
Although the constraints from this initial Rocky Planet Search observing run are relatively mod- 



^Detectability has been extrapolated to the masses in this range higher than were simulated, however this contir- 
butes a negUgible uncertainty as the detectabilities for these more massive mass planets are very close to 1.0 
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est, our simulations offer a direct metfiodology to determine an empirical exoplanet mass function. 
They can readily be improved with (1) an expanded sample of high precision data on a larger 
sample, (2) simulations covering larger sets of high-precision AAPS matching the quality obtained 
in this intensive observing run, (3) higher resolution simulations (i.e. more realisations at more 
closely sampled orbital parameters) and (4) further improvements in our understanding of stellar 
jitter and convective granulation. 
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prove the paper. We acknowledge support from the following grants: NSF AST-9988087, NASA 
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Fig. 1. — Histograms of measured Dopplcr velocities for eaeh of our 24 targets. Also plotted are 
Gaussians with FWHM equal to the stellar jitter (dotted) and the velocity scatter (dashed). His- 
tograms of the residuals are shown for the two objects with known planets (HD 4308 and HD 16417), 
while the binary system a Cen (HD 128620 and HD 128621) have had a linear trend subtracted. 
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Fig. 2. — Data from Rocky Planet Search run for HD 84117. 
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Fig. 3. — 2DKLS Periodogram evaluated at e = 0.44 determined for HD4308 from AAT Rocky 
Planet Search data alone. 




Fig. 4. — Keplerian fit to HD4308 data from Rocky Planet Search (RPS) run. 
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Fig. 5.— 2DKLS periodogram determined at e = 0.24 for HP 4308 from AAT da ta (2005 Oct-2007 
Nov, including Rocky Planet Search data) and HARPS data (jUdrv et al.ll2006l l. The inset shows 
2DKLS power as a function of eccentricity at constant period, and demonstrates how weakly this 
data set (in common with most Doppler data sets) constrains eccentricity. 
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Fig. 6.— Velocities for HP 4308 from the AAT (2005 Oct-2007 Nov, including Rocky Planet Search 
data) and HARPS (jUdrv et al.ll2006l ) plotted against time (upper panel) and phased at the best-fit 
orbital period (lower panel). The line indicates the best-fit Keplerian to this combined data set, 
the parameters of which are listed in Table |H 
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Fig. 7.— 2DKLS at e=0.0 for the HD 4308 velocities of Table [3] after removal of the best fit 
Keplerian of Fig [6l While there is some evidence for excess power at 30-80 d periods, this is not 
currently sufficient to justify claiming the detection of a second planet. (The vertical dashed line 
indicates the time span of the observations used to make this periodogram.) 
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Fig. 8. — Detectability contours for planets as a function of input period and input planet mass (i.e. 
Doppler m sin i minimum mass) for two example stars (HD 4308 and HD 53705) , and averaged over 
all 24 stars. (See text for our formal definition of "detectability"). The known exoplanet HD 4308b 
is shown as an asterisk. The data for HD 4308 allow us to recover the previously detected planet - 
in comparison the HD 53705 data demonstrate significantly larger sensitivity to lower-mass planets, 
due to the better sampling that was achieved on this run for this star. Averaged over the whole 
sample, the 10% detectability contour extends down to 3 for periods from 2-10 d. 
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Fig. 9. Expected number of exoplanet detections in our survey of 24 stars for different mass 
functions. The asterisks with uncertainties indicate the exoplanets HD 16417b and HD 4308b with 
msmi masses of 0.0692 and 0.0442 Mjup respectively. 
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Not e. — Coordinates and magnitudes a r e based on the Hipp a rcos catalogue ( Ferryman et al.M 19971 ): spectral types are from iHouk 
1 19821 ): activities are from lHenrv et all (|l996l ): pTinnev et al.l l|2002b ): Ijenkins et al.l (|2006l ): metallicities. isochrone masses and ages are from 



Valenti fc Fischej l|2005l . VF05) and iNordstrom et all (|2004l . N04). 
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Table 2. Rocky Planet Search Observations 
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Note. — A^obs indicates the number of nights on which 
the object was observed. An observation on any particular 
night always comprises at least three separate exposures. Also 
given is the velocity scatter (RMS), median velocity uncer- 
tainty, activity jitter (see tex t), and oscillation jitter (from 
O'Toole. Tinnev fc Jones 20081 ) . The latter is based on the av- 



erage exposure time for each target and is given as a guide. 
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Table 3. AAT Doppler velocities for HD 4308. 



JD Velocity JD Velocity 

(-2450000) (ms-^) (-2450000) (ms"^) 



3665.119898 


-0.3 


± 


0.9 


3666.110456 


-0.6 


± 


0.9 


3668.136398 


2.8 


± 


0.8 


3669.119458 


0.6 


± 


0.9 


3670.126654 


4.5 


± 


1.1 


3671.147771 


3.6 


d= 


1.3 


3698.032791 


-9.0 


± 


1.0 


3700.063296 


-0.7 


± 


0.8 


3702.062004 


1.1 


± 


1.3 


3938.327905 


5.9 


± 


1.0 


3942.284604 


-0.1 


± 


0.9 


3943.328172 


2.6 


± 


0.8 


3944.283860 


-1.3 


± 


0.9 


3945.320288 


-3.8 


± 


0.9 


3946.318428 


-4.6 


± 


0.9 


3947.328379 


-5.5 


± 


1.0 


4008.152191 


-1.1 


± 


1.0 


4009.177898 


-0.0 


± 


1.0 


4010.168991 


0.2 


± 


0.9 


4011.141538 


-2.4 


± 


0.9 


4012.129794 


-2.5 


± 


0.9 


4013.171260 


-0.2 


± 


0.9 


4014.188894 


-2.6 


± 


1.2 


4016.210211 


-1.5 


± 


1.0 


4018.142500 


-0.1 


± 


1.0 


4037.046606 


8.5 


± 


0.8 


4038.125336 


4.6 


± 


1.0 


4039.054689 


0.4 


± 


1.0 


4040.053384 


-2.3 


± 


0.9 


4041.068042 


-0.4 


± 


0.8 


4111.941520 


3.1 


± 


1.2 


4112.942694 


4.1 


± 


1.1 


4113.949634 


0.3 


d= 


1.7 


4114.947273 


-0.4 


± 


1.0 


4118.926666 


-8.9 


± 


1.1 


4119.920872 


-6.7 


± 


1.1 



4120.935790 


-2.4 


± 


1.7 


4125.954730 


1.7 


± 


1.1 


4126.924257 


2.6 


± 


1.2 


4127.931946 


3.7 


± 


1.1 


4128.927984 


0.5 


± 


1.0 


4129.930445 


1.1 


± 


1.0 


4130.924765 


-0.1 


± 


0.9 


4131.937199 


-0.2 


± 


1.1 


4133.927357 


-2.3 


± 


0.9 


4134.931174 


0.5 


± 


1.0 


4135.948345 


-1.5 


± 


0.9 


4136.945095 


3.1 


± 


1.0 


4137.977861 


6.6 


± 


1.3 


4140.931894 


7.5 


± 


1.3 


4145.918589 


5.0 


± 


1.2 


4149.900027 


-4.6 


± 


1.1 


4153.907836 


1.6 


± 


1.1 


4154.918636 


4.5 


± 


1.2 


4253.327811 


0.4 


± 


0.9 


4255.246671 


-3.1 


± 


1.3 


4256.321705 


-2.6 


± 


1.2 


4257.304964 


-3.6 


± 


1.7 


4334.157294 


3.0 


± 


1.2 


4335.254660 


7.8 


± 


2.7 


4336.264333 


-3.7 


± 


1.3 


4338.306202 


-8.4 


± 


1.4 


4369.133832 


-0.3 


± 


0.9 


4371.241454 


-1.5 


± 


1.0 


4372.138956 


5.5 


± 


1.2 


4373.185170 


2.7 


± 


1.3 


4375.180751 


0.6 


± 


1.0 


4425.078066 


0.2 


± 


1.1 


4429.012543 


-2.0 


± 


1.1 


4429.946086 


-4.9 


± 


0.9 


4431.004964 


-6.9 


± 


1.2 


4431.997618 


-5.8 


± 


1.2 



Note. — Julian Dates (JD) are heliocentric. Velocities 
are barycentric but have an arbitrary zero-point. 
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Table 4. Orbital parameters for HD 4308b. 



Parameter 


AAT & HARPS 


HARPS 


Orbital period P (days) 


15.609 ± 0.007 


15.56 ± 0.02 


Semi-amplitude K (ms~^) 


3.6 ± 0.3 


4.07 ± 0.2 


Eccentricity e 


0.27 ± 0.12 


0.00 ± 0.01 


Periastron (JD- 2450000) 


108.5 ± 1.9 


3314.7 ± 2.0 




210 ± 21 


359 ± 47 


msini (M^) 


13.0 ± 1.4 


15.1 ± 0.5^^ 


Semi-major axis (AU) 


0.118 ± 0.009 


0.119 ± 0.009*^ 


Nfit 


113 


41 


RMS (ms-i) 


2.6 


1.3 




1.17 


1.3 



^msini and semi-major axis numbers are those derived from 
the kinematic parameters of lUdry et al.l (120061 ). but assuming the 
same host-star mass (0.91±0.05 Mq) used in this paper. 



